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the Brillouin improved Wilson action, we are able to calculate this ratio in a single framework, 
using a relativistic fermionic action throughout. The calculation is carried out on selected ensem- 
bles of two flavor clover improved lattices produced by the QCDSF collaboration, allowing an 
extrapolation to the continuum, to infinite volume and to the physical pion mass. 
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1. Introduction 

In this presentation we detail our measurement of the charm to strange quark mass ratio 
nic/ms [1]. The exact determination of the masses of quarks is of particular phenomenological 
interest, since they cannot be measured directly in experiments, leading to rather large errors in 
their PDG values [2]. This is an area where lattice QCD can provide important input, especially 
in the case of ratios of quark masses where lattice renormalization factors cancel, which otherwise 
require dedicated renormalization techniques to keep the systematic error small. 

We use the Brillouin action, which we introduced in Ref. [3], to calculate mc/m,. Among a 
number of nice features, the Brillouin action has been shown to exhibit small cut-off effects even 
in the charm quark mass region when compared to regular clover-improved Wilson. In the sea 
sector, we select a set of 2-flavor clover-improved lattice gauge ensembles provided by the QCDSF 
collaboration which allow us to extrapolate to zero lattice spacing, to the physical pion mass and 
to infinite volume. We compare our result, and find agreement to the two recent state-of-the-ait 
calculations of this quantity by the HPQCD collaboration [4] and the ETM collaboration [5]. 



2. Lattice Setup 

In this work we use a mixed action approach to determine the ratio mc/m^. In the sea sector, 
we use configurations generated by the QCDSF collaboration. The parameters of the gauge con- 
figuration ensembles used ai^e detailed in Table 1, while for the simulation details see Refs. [6, 7]. 
The selection of the specific subset of configurations was chosen such as to provide a wide enough 
window in the pion mass, the physical volume and the lattice spacing in order to allow for an ex- 
trapolation to the physical point. In the valence sector of our mixed action approach we use our 
recently developed Brillouin improved Wilson action [3]. As in our original publication of this 
fermion action, we perform a single step of APE-smearing on the gauge links to further improve 
gauge-noise and the approach to the continuum, as well as tree-level clover improvement {csw = 1)- 

As will be more apparent in the following sections, in our approach we need only set the scale 
in the final step, when performing the extrapolation to the physical point. In that step, the scales 
used are taken from Ref. [7], using the Sommer scale as ro = 0.5 fm at the chiral limit. This yields 
a = 0.076, 0.072 and 0.060 fm for (5 = 5.25, 5.29 and 5.4 respectively. 



3. Methodology 

Our mixed action approach requires us to tune the Wilson hopping parameter for the strange 
and charm quarks. For this tuning we choose ratios of meson masses, which allows tuning these 
quantities without the need to set the scale. Namely, we choose the two ratios: 

7?, = ^%^ and R, = ^-^^. (3.1) 

Note that both ratios depend on both the strange and the charm quark mass. The procedure we 
follow is to invert for a few choices of the two hopping pai^ameters and Kc, and to compute R\ 
and /?2 for every combination of these two. Typically we invert for three and three Kc values, 
which leaves us with nine points. We then fit to find two curves in the (k;, fQ.) coordinate system 
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Table 1: Table of the QCDSF configurations used. For each ensemble ^500 configurations were used. 
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Figure 1: Example of tuning for and (Q . We invert for three values of fc, and fQ giving us nine R\ (left) 
and /?2 (right) ratios (blue circles). The blue surfaces are fits to these nine points. The red surface shows the 
constant surface of the physical value of the respective ratio, while the red curve is the intersection of the red 
and blue surfaces (shown as a dashed line on the base-surface of each plot). The point where the two curves 
intersect is taken as the tuned K"-pair for this ensemble. 



along which R\ and R2 acquire their physical values. We take these as R\ = 0.801 and R2 = 12.402. 
This gives us a single intersection point, and the K"-values at this point we call the "tuned" kappa- 
values: K* and K*. An example of this tuning procedure for one of the ensembles considered is 
shown in the plots in Fig. 1 . 

For this calculation, we use the Partially Conserved Axial Current (PCAC) relation: 

where /, 7 = c or s. The ratio of quark masses is then given by the ratio of Axial Ward Identity 
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quark masses which are associated with the PCAC quark masses through: 

AW/ ^PCAC 

= [1 +ci{mY- mJ){bA - bp) + ff{a^)] . (3.3) 

nij rrij 

At tree-level, we have ca = and bA — bp = 0, meaning we can take directly the quark mass ratio 
r = mc/nis from the ratio of PCAC quark masses. A more detailed discussion of these terms is 
given in Ref. [1]. Additionally, we avoid involving correlators of two charm propagators by using 
the relation: 

Ms 

This way the mf^"^*- mass never appeai^s in our analysis. To summarize, the procedure we follow is 
to first determine the "tuned" K"-pairs as we have akeady detailed above. Subsequently, we compute 
the ratio r for each of the nine (fc,, Kc) pairs, as in the case of /?i and R2, and by fitting to a surface 
we can interpolate for the value of the ratio r at the tuned point ( K"* , K"* ) . This is taken as the value 
of nic/nis for the given ensemble. 




Figure 2: Example of the determination of 
mc/rrig at the tuned K values. The blue surface 
is a fit to the blue circles. The dashed lines cross 
at the tuned K values, as shown in Fig. 1. The 
black asterisk is the central value of nic/ms for 
this ensemble. This is done for each Jack-Knife 
sample to obtain a Jack-Knife error for the ratio. 
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Figure 3: Example of the error associated with 
the determination of the tuned K values. The plot 
is a top-down view of Fig. 2, with error bands 
added to the dashed lines. The error bands where 
obtained through a Jack-Knife analysis. 



This entire procedure is carried out within a Jack-Knife analysis in order to propagate the 
statistical error to the final value of nic/nis. In Fig. 2 we show an example of the interpolation to 
the tuned K"-values. In Fig. 3, we show the statistical errors associated in determining the tuned K- 
values. These are obtained by performing a Jack-Knife analysis, where the crossing is determined 
in every Jack-Knifed sample. We note that the errors on K* and (C* are not needed in the final 
analysis, since the Jack-Knife procedure is carried out up to the final determination of r. The eiTors 
in Fig. 3 have been extracted for illustration purposes only. 
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Figure 4: Examples of fits of our data to the global fit function Eq. 4. 
independent variables, in each of these three figures, the data points are i 
of the two axes not shown. The left plot shows a fit to f^''^'^{a,Mjc,L) as a function of a^, the center plot 
a fit to f^-^'^{a,Mj,:,L) as a function of the pion mass squared and the right plot a fit to f^'^-^ {a,M:^,L) as a 
function of the inverse volume. 



4. Results 



As mentioned, the ensembles were chosen in order to allow for a reliable extrapolation to 
the physical pion mass, infinite volume and zero lattice spacing. We perform a global fit so as to 
extrapolate simultaneously to these three physical limits. For each of the three extrapolations we 
take two fit ansatze which means these can be combined in eight ways, thus yielding eight global 
fit functions: 



f''J'\a,M„,L) = Co \+c\F'{a) +c{&{M„)+c\h\M„,L 



with /,;,/t = Oorl (4.1) 



where: 



F^\a) 



■Ml 



F\a) 
G\M„) 

M\M„,L) 



-Ml 



(4.2) 



In Fig. 4 we show three plots of representative extrapolations. In each plot, the extrapolation is 
shown with respect to one of the three extrapolation axes. For illustration purposes, the data in the 
plots have already been extrapolated for the axes not shown. 

We note here that we considered two possible ways to discretize the partial derivative involved 
in the calculation of the PCAC quark mass (Eq. 3.2). Namely, we considered the regular symmetric 
derivative: dtf{t) = [/(? + 1) — f{t — l)]/2 as well as the improved form: 



fit - 2) - 8/(? - 1) + 8/(? + 1) - /(? + 2) 



12 



(4.3) 



Although the two forms yield consistent results for mc/ms, using the latter form we obtain ;^^/d.o.f 
of order one for each of the eight extrapolations, while using the regular form we obtain ;^^/d.o.f in 
tendency above one. We therefore opt to use the improved form as the definition of the discretized 
derivative dt throughout our analysis. 
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Figure 5: The ratio mc/nis as obtained from the 
eight fit functions. The solid vertical line is the 
average of the eight values. The lighter gray 
band indicates the statistical error of the central 
value while the darker gray band shows the sys- 
tematic error, obtained from the standard devia- 
tion of the eight values. 



Figure 6: The "control ratio" Rt, obtained by 
the eight fits, normalized by its physical value 
of 1.7707. 



The fit procedure described gives eigiit values of tlie ratio md nis at zero lattice spacing, infinite 
volume and the physical pion mass, one for each of the eight [i,j,k] combinations of Eq. 4.1. We 
combine these eight values to obtain our final estimate of the ratio of mc/m^ , by taking the standard 
deviation of these eight values as the systematic eiTor of our measurement, while the average of 
the eight central values and of their statistical error is taken as the central value and statistical error 
respectively. The eight results of the fit ai^e shown in Fig. 5 with the error-bands denoting the 
systematic and the statistical error. Our final value for the ratio of rtic/ms is: 

m,/m, = 11.34(40)stat(21)syst. (4.4) 

Ml 

As a potential check of our method, we carry out an identical analysis for the ratio = y 

of which the physical value is well determined at 1.7707. The result of the eight fits for this 
"control quantity" is shown in Fig. 6 (normalized to the physical value). The result is = 
1.77(ll)stat(07)syst- The excellent agreement of this quantity with its physical value indicates that 
our tuning procedure and analysis are sound and that the obtained estimate of mc/nis is reliable. 



5. Summary 

In this presentation we have carried out a determination of the fundamentally interesting quan- 
tity of the charm to strange quark mass ratio mc/nis. This was done using a relativistic fermion 
action in the valence sector with improved scaling at large quark masses. We perform a controlled 
extrapolation to the physical pion mass, zero lattice spacing and infinite volume. 

Our result is in reasonable agreement with two other recent calculations of this quantity. 
Namely, the first result obtained using a relativistic fermion action was a calculation using stag- 
gered quarks [4] which aiTived at mdms = 11.85(16), while a second result available by the ETM 
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HPQCD (HISQ) 
H-H ETMC (N/=2 twisted mass) 
^ This work 

11.0 11.5 12.0 12.5 

Figure 7: Three results of the ratio mc/nis computed on the lattice. The red star corresponds to the result by 
the HPQCD collaboration, the blue square is the result by the ETM collaboration and the green circle is the 
result of this work. 

collaboration [5] arrived at mc/nis = 12.0(3), wliere tlie latter's error is statistical witli no clear 
indication of tlie systematic uncertainty involved. A comparison of our result to these two de- 
terminations is shown in Fig. 7. Although our result is somewhat more imprecise, it serves as 
an important check since our fermion action is free from any isospin or taste symmetry breaking 
which is inherent in either staggered or twisted mass fermions. 

Additionally, we note that the computer resources required for this calculation where rather 
modest. The improved approach to the continuum exhibited by the Brillouin improved fermion 
action allowed a reliable extrapolation to the physical point with a moderate number of statistics 
and gauge ensembles. Namely, the requirements for tiiis calculation where roughly 20 node-years, 
of an eight-core Nehalem node. This number includes inverting for three kappa values for each of 
the charm and strange quark masses, as well as trial inversions at each ensemble to determine the 
best three such values. 
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